نشریه مهندسی عمران فردوسی سال بیست و پنجم شماره یک, ۱۳۹۲ 


به کار گیری دبنامیک سیالات محاسباتی در مدل‌سازی سلول‌های جربان انویه در کانال ذوزنقه‌ای * 


۱ مهنا تاج‌نسایی"۳ تفا گیوهخی 


مهدی اژدری مقدم!" 
چکیده جریان عبوری ا زکانال‌ها دارای سه مولفه‌ی سرعت, یک مولفه در جهت جریان و دو مولفه در جهت عرض ی کانال است. در اثر 
ناهمکنی نوسانات سرعت, یک سر ی گردابه‌های چرخحشی در مقط عکانال ایجاد می‌گردند که سلول‌های جریان ثانویه امیده می‌شوند. این 
سلول‌ها سبب ایجاد تنش برشی عرضی در جداره‌های کانال می‌شوند و در نتیجه بحث فرسایش وآبشستگی جداره‌های کانال را مطرح 
می‌کنند. ناهمگن یآشفتگی به‌علت شرایط مرزی بستر دیواره‌ی جانبی و سط حآزاد, نسبت ابعادی کانال و هندسه‌ی کانال به‌وجود م یآید. 
مکانیزم ایجاد جریانهای ثانویه را می‌توان به‌کمک معادله‌ی چرخحش طولی بیان نمود. در اين تحفیق با استفاده از دینامیک سیلات محاسباتی 
(رل) و نرم‌افزار ) ۸515 جریان در کانال‌های ذوزنقه‌ای مدل‌سازی خحواهد شد و علاوه بر بررسی عملکرد 2 آن) در مدل‌سازی 
سلول‌های جریان ثانویه, توزیع سرعت عمق میانگین و تنش برشی مرزی د رکانال مورد بررسی قرار خواهد گرفت. در این راستا برای 
ارزیابی مدل ساحته شده از مدل آزمایشگاهی تومیناگا و همکاران و مدل تحلیلی شینو و نایت بهره گرفته شده‌است. با صحت‌سنجی مدل 
عددی ساخته‌شده براساس نتایج آزمایشگاهی توافق خوبی مشاهده گردید. نتایج حاکی از عملکرد به‌تر مدل آشفتگی تنش برشی رینولدز 5560 

در مدل‌سازی سلول‌های جریان بود. سطح همگرایی ۱۰ به‌دلیل پایداری نتایج در آن به‌عنوان سطح همگرایی بهینه انتخاب شد. 


واژه‌های کلیدی کانال ذوزنقه‌ای. سلول‌های جریان‌انویه. دینامیک‌سیالات‌محاسباتی ((۳1:» سرعت عمق‌میانگین. توزیع تنش‌برشی- 


مرزی. 
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,((۳۲1) قتصجمصنا فتبطط امممتام‌و‌وم روااعی مت وتهل‌همععو. رقاعصصهطن اهلزم2ع۲2۵ ۲ 16۷۷۵۲۵۵ 
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* تاریخ دریافت مقاله ۸۹/۸/۲۵ و تاریخ پذیرش آن ٩۱/۱۱/۱‏ می‌باشد. 
(۱) نویسنده مسژول دانشیان دانشگاه سیستان وبلوجستان, دانشکده مهندسی, گروه مهندسی عمران. 
(۲) عضو هیات علمی دانشگاه آزاد اسلامی واحد بیجار بیجار. ایران. 


(۳) استادیا دانشگاه سیستان وبلوجستان, دانشکده مهندسی گروه مهندسی عمران. 


‌ 


مقدمه 

با توجه به اهمیت کانال‌ها در بحث انتقال آب و 
وابستگی اعتبار سیستم‌های مدیریت آبی به میزان توجه 
آن‌ها به انتخاب شرایط کانال‌های انتقال آب. بررسی 
شرایط جریان در آن‌ها از اهمیت ویژه‌ای برخوردار 
است» به‌گونه‌ای که کاهش ظرفیت انتقال آب به کاهش 
منابع آب قابل دسترس و در نتیجه افت قابل توجه 
منابع اقتصادی منجر می‌شود. از مهم‌ترین موضوعات 
بررسی جریان در کانال‌هاء می‌توان به اثر مولفه‌های 
عرضی سرعت روی جریان اشاره نمود. که بارزترین 
این اثرات ایجاد سلول‌های جریان انوبه است. در حال 
تحار روفن‌های که بزای شعاشیه وطر این کانال‌هتا 
وجود دارد اثر جریان‌های ثانویه و در نتیجه توزیع 
نیروهای هیدرودینامیکی موجود روی پیرامون مرطوب 
را نادیده می‌گیرند و این مسأله یکی از مهم‌ترین دلایل 
تغییر شکل مقطع عبوری» کاهش قابلیت اعتماد 
عمل کرد و ظرفیت کانال‌ها است. 

در این زمینه مطالعات متعددی به‌صورت 
آزمایشگاهی و عددی صورت پذیرفته‌است که در این 
میان به‌صورت نمونه می‌توان به مواردی اشاره نمود. از 
جمله: بررسی توزیع تنش برشی مرزی در کانال‌های 
منشوری و داکت‌ها به‌صورت تجربی توسط نایت و 
دمتریو [۱) نایت و همکاران [2] لای و نایت» رودس 
و نایت [۱] و یانگ و لیم [3]» تومیناگا و همکاران[4]» 
جابرزاده و همکاران [۵] ابداع روش شینو و نات 
7 1۷۲۵۵۰00 اطم‌نیکا 0ج 0ط60) برای محاسبه‌ی 
توزیع‌های جانبی سرعت و تنش برشی عمق میانگین 
برای جریان در کانال‌های منشوری مستقيم کاربرد 
روش 5101 برای کانال‌های مختلف توسط نایت و 
همکاران [6) خداشناس و همکاران [7]براق جزفیات 
بیش‌تر از منبع [۱] می‌توان بهره گرفت. 

از آن‌جایی که کانال‌های با مقطع ذوزنقه‌ای 
معمول‌ترین شکل مقطع برای انتقال آب به شمار 


می‌روند. در این مقاله به بررسی سلول‌های جریان 


نشریه مهندسی عمرآن فردوسی 


بهکا رگیری دینامیک سیالات محاسباتی در مدل‌سازی ... 


ثانویه در کانال‌های مستقیم با مقطع ذوزنقه‌ای پرداخته 
شده است. 

مانند ساير پدیده‌های فیزیکیء برای بررسی 
سلول‌های جریان انویه نیز می‌توان از هرسه نوع 
دینامیک سیالات تجربی. تثوری و محاسباتی بهره 
جست. اطلاعات دقیق در مورد یک فرآیند فیزیکی 
غالباً توسط اندازه‌گیری عملی به‌دست می‌آید. یک 
پیش‌گویی تشوری حداکثر استفاده را از نتایج مدل 
ریاضی خواهد کرد و در مقایسه با آن نتایج تجربی را 
مورد استفاده کم‌تری قرار می‌دهد. در این زمینه 
فعالیت‌های گسترده‌ای در قالب مطالعات آزمایش گاهی 
و به‌ویزه تحلیلی صورت گرفته است. پیچیدگی 
معادلات حاکم بر مسأله تأثیر متقابل فیزیکی مختلف؛ 
گذرا بودن اغلب مسائل مهندسی. بالا بودن هزینه‌های 
مربوط به تجهیزات آزمایشگاهی و محدودیت استفاده 
از دستگاه‌های اندازه‌گیری در بسیاری از کاربردهای 
عملی: ازخمله بدلایلن. هسد که استفاده آودروش‌هناق 
تحلیلی و آزمایشگاهی را در مقایسه با روش‌های 
عددی محدود می‌کند [۸]. 

در این مقاله به‌ کمک دینامیک سیالات محاسباتی 
(0۳0) و با استفاده از نرم‌افزار محاسباتی ۸5۷5 
نسخه‌ی 12.0 شکل‌کیری سلول‌های جریان 
انویه در کانال‌های ذوزنقه‌ای مورد بررسی قرار گرفته 
و برای ارزیابی مدل ساخته شده در محیط نرم‌افزار از 
مدل آزمایشگاهی تومیناگا و همکاران [4] کمک گرفته 


شده است. 


سلول‌های جریان ثانوبه 
جریان‌های انویه جریان‌هایی هستند که در صفحه‌ی 
مقطع جریان و يا حول محوری عمود بر صفحه‌ی 
مقطع جریان به‌وجود می‌آیند. این گونه جریان‌ها در 
انحناها به جریان‌های ثانوبه‌ی قوی موسوم می‌باشند 
ولی در کانال‌های بدون انحنا در پلان نیز به‌جهت تأثیر 


سال بیست و پنجم شماره یک ۱۳۹۳ 


مهدی اژدری مفدم- مهنا ناج نسایی- محم دکیوه‌چی 


زبری جدار و نامنظمی مقطع. جریان‌های انویه‌ی 
ضعیف ایجاد می‌گردند .]٩[‏ 

به طورکلی می‌توان جریانات چرخشی رابه سه 
نوع اصلی چرخش حول محور عمقی محور عرضی و 
محور طولی تقسیم‌بندی کرد [۵]. جریان‌های انویه 
نتیجه‌ای از ناهمگنی آشفتگی هستند که به‌دلیل شرایط 
مرزی بستر دیواره‌ی جانبی و سطح آزاده نسبت 
ابعادی کانال و هندسه‌ی کانال ایجاد می‌شوند. این نوع 
جریان‌ه | نسبت به هندسه‌ی کانال بسیار حساس 
هستند. تومیناگا و همکاران بررسی جریان‌های ثانوبه 
در مهندسی هیدرولیک در جریان کاناهای روباز را 
به‌دلیل اثر این جریان‌ها روی میدان جریان میانگین 
اصلی. قانون اصطکاک. تشکیل اشکال سه‌بعدی بستر 
از قبیل نوارهای شنی و انتقال رسوب سه‌بعدی با 


اهمیت می‌دانند [4], 


مکانیسم تولید جریان‌های ثانوبه 

علت ایجاد جریان‌های عرضی در جریان طولی 
یکنواخت. ناهمسانی تنش‌های رینولدز و به‌طور 
اساسی مژلفه‌های عمودی آن است. میدان سرعت 
میانگین طولی به‌وسیله‌ی تنش‌های برشی آشفته تشک 
می‌شود. اما جریان‌های عرضی اآثر قابل ملاحظه‌ای 

پرانتل دو مکانیسم اصلی به‌وجود آورنده‌ی 
جریان‌های چرخشی در جهت جریان را به‌شکل زیر 
بیان کرد [۱۱]: 
۱- انحراف تنش‌های برشی اصلی به‌وسیله‌ی یک 
گرادیان فشار عرضی و یانیروی جسمی (مانند 
کانال‌های منحنی شکل) 
۲- ناهمسانی و نوسانات 1 شفک در صفحه‌ی عمود بر 
محور جرخحش. 
مکانیسم تولید جریان‌های انویه‌ی القا شده توسط 
آشفتگی در جریان کانال مستقيم به‌کمک معادله‌ی 


سال بیست و پنجم. شماره یک ۱۳۹۲ 


چرخش طولی توضیح داده می‌شود. معادله‌ی چرخحش 
طولی در جریان آشفته کاملاً توسعه یافته به‌صورت زیر 


داده می‌شود 41 
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که در آن ۷۷ و ۷ به‌ترتیب مولفه‌های سرعت در 
تور ای ی نم سس کافال ۳ ۳) وا 
به‌عبارتی مژلفه‌های سرعت جریان ثانویه هستند و 


به‌کمک تابع پتانسیل ۷ به‌صورت زیر بیان می‌شوند: 


۳ كٍِِ-۷ , ۷۷ 
0۷ 02 


که در آن 9 ویسکوزیته‌ی سینماتیکی سیال» 
تی ری و ۱ متام 
عمودی رینولدز می‌باشند. توزیع تفاضل تنش‌های 
عمودی ( 2 قب) نقش اساسی و مهمی را در ساختار 
جریان‌های ثانویه بازی می‌کند. به گونه‌ای که تفاوت 
بین الگوی جریان‌های انویه در کانال‌های روباز و 
بسته به کمک توزیع تفاضل تنش‌های عمودی 
( 2 قب) توضیح داده می‌شود. در شکل (۱) توزیع 
این تفاضل در سه نوع کانال نمایش داده شده است. 
جرارد نیز این توزیع را تعیین کننده‌ی ساختار 
جریان‌های انویه می‌داند [4]. 


نشریه مهندسی عمرآن فردوسی 
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بهکا رگیری دینامیک سیالات محاسباتی در مدل‌سازی ... 


1 برد 68 9 


شکل ۱- خطوط بی‌بعد پبر۲2(/272- 0۳2 9 پر۲۳2(/72 - 72) از راست 0)کانال بسته ۲) کانال ذوزنقه‌ای ۳)کانال روباز[ ی ] 


اثرات سلول‌های جریان ثانویه روی جریان 
اصلی 

سرعت ماکزیمم در هر مقطع قانم. در توا نکب سطح 
آب و در فاصله‌ی ۵ ] ۰۰*۰۵ عمق جریان از سطح 
آزاد اتفاق می‌افتد. علت اصلی این که سرعت ماکزیمم 
در نزدیکی سطح آزاد پیش می‌آید بیش از آن‌که تصت 
تأثیر تنش برشی ناشی از مقاومت هوا باشد» تحت تأثیر 
جریان‌های ثانویه ضعیف می‌باشد ]٩4[‏ 

جریان‌های عرضی (ثانویه) ب‌صورت قابل توجهی 
فرآیندهای انتقال احتلاط و گرما و جرم را افزایش می- 
دهند. نیروهای هیدرودینامیکی را روی پیرامون 
مرطوب کانال بازتوزیع می‌کنند» نواحی با کسر 
رسوبات ایجاد می‌کنند و به موجب آن وقوع انبساط 

وجود جریان‌های چرخشی در مسیر یک رودخانه 
رودخانه می‌شود [۵]. جریان‌های انویه. ساختارهای 
سه‌بعدی ایجاد جریان میانگین اصلی را تحت تأثیر قرار 
می‌دهد [4]. نتایج بررسی‌ها نشان دادند که سلول‌های 
جریان انویه روی توزیع تنش برشی مرزی تأثیر 
می‌گذارند [12]. 


دینامیک سیالات محاسباتی 


نشریه مهندسی عمرآن فردوسی 


پدیده‌های همراه نظیر واکنش‌های شیمیایی براساس 
شبیه‌سازی کامپیوتری است [۸ و ۱۳] پیچیدگی 
معادلات حاکم بر مسأله‌ی تأثیر متقابل فیزیکی مختلف؛ 
گذرا بودن اغلب مسائل مهندسی, بالا بودن هزینه‌های 
مربوط به تجهیزات آزمایشگاهی و محدودیت استفاده 
از دستگاه‌های اندازه‌گیری در بسیاری از کاربردهای 
عملی: از جمله دلایلی خستد که استفاده از زوش‌های 
تحلیلی و آزمایشگاهی را در مقایسه با روش‌های 
عددی محدود می‌کند [۸. 

در نهایت یک برنامه‌ی ۳10) خوب شامل تعادل 
بین حل دقیق و پایداری جواب می‌باشد. بررسی اعتبار 
10 نیازمند اطلاعات بالا در جزئیات مربوط به 
شرایط مرزی مسأله و به‌دست آوردن حجم بالایی از 
نتایج است. برای بررسی اعتبار از یک راه معنی‌دار 
ضروری است که داده‌های تجربی در یک میدان تشابه 


تاش ان ۱۲۱ ] 


مدل‌سازی عددی سلول‌های جریان ثانوبه 
اندازه گیری مستقیم جریان‌های انوبه در کانال‌های 
روبان با توجه به اين نکته که سرعت جریان‌های ثانوبه 
درصد کمی از سرعت طولی جریان است. کار دشواری 
است [4] و نیاز به تجهیزات آزمایشگاهی دقیق و دقت 
بالای کاربر دارد. اگرچه محققان تلاش‌های بسیاری 
برای شبیه‌سازی تنش برشی مرزی بدون مراجعه به 


0 داشته‌اند. آن‌ها خیلی موفق نبوده‌اند [14]. برای 
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محاسبه‌ی اثر جریان‌های عرضی روی فرآیندهای 
مختلف. نخست نیاز است که یک مدل ریاضی از 
میدان سرعت میانگین سه‌بعدی در شکل‌های مختلف 
مقطع عبوری جریان شامل مرزهای ناهمگن با توانایی 
مدل‌سازی فرآیندهای جریان در بازه‌ی گسترده‌ای از 
ناحیه‌ی تغییر محبط‌های کانال ایجاد کنیم. در مرحله‌ی 
دوم بررسی‌های آزمایشگاهی و تجربی در محل انجام 
شود مشخصه‌های اصلی توزیع خصوصیات آشفتگی 
سه‌بعدی معلوم شود و داده‌های تجربی برای مستند 
ساختن مدل‌های ریاضی ایجاد شده به دست آید [10]. 

برای بهره‌گیری از ۳10 و مدل‌سازی عددی. نیاز 
به استفاده از نرم‌افزارهایی برای رسم هندسه 
شبکه‌بندی» معرفی شرایط جریان, تحلیل جریان و 
مشاهده‌ی نتایج وجود دارد. برای رسیدن به این اهداف 
دو نرم‌افزار 0۳5 10 و 0۳ ۸5۷5 انتخاب 
شده‌اند. در ادامه روند مدل‌سازی به‌ترتیب ارائه 


می‌شود. 


آشنایی کلی با نرم‌افزار 0۳۲ ۸۷55. نرم‌افزار 
۶ ۵5۷5 یک نرم‌افزار با هدف عمومی دینامیک 
سیالات محاسباتی است که یک حلگر پیشرفته را با 
قابلیت‌های پیش و پس‌پردازشگر قدرتمند ترکیب 
کرده‌است و توانایی مدل کردن مواردی از قبیل 
جریان‌های دائمی و غیردائمی» جریان‌های آرام و 
آشفته, جریان‌های مادون صوت و در حد صوت و 
مافوق صوت. نیروی شناوری. جریان‌های غیرنیوتنی؛ 
انتقال مولفه‌های اسکالر فاقد عکس‌العمل. جربان‌های 


چندفازی» مسائل احتراق و ... را دارا می‌باشد [12]. در 
نرم‌افزار 0۳26 ۵5۷5 چهار سطح همگرایی مقادیر 
کته از 10۳ ۳ 10 
مقادیر کم‌تر از *107 تعریف شده است. 

در اين نرم‌افزار چهارده مدل برای مدل‌سازی 
آشفتگی وجود دارد شامل مدل ع- ۰6 مدل انتقال 
تنش برشی (1120500۲1 ووع50۲ 5۳68۲) با 991 . مدل 
تنش رینولدز ۹50 (50655 6700105 550). مدل 
۲-0 و ... می‌باشد. نرم‌افزار 0۳ ۸۱۷5۷5 روش 


عددی حجم محدود را مورد استفاده قرار داده و ز‌ 
الگوریتم 00100160 برای گسسته‌سازی معادلات بهره 


می گیرد. 


مدل‌سازی هندسی کانال. برای مدل‌سازی هندسه‌ی 
کانال موردنظر دو نرم‌افزار 028 ماند۸ و 0۳۳ 10516 
مورد ارزیابی قرار گرفته و از بین این دو نرم‌افزار 
گزینه‌ی برتر نرم‌افزار 0۳ 1015۷۲ به‌دلیل توانایی آن 
هم در رسم هندسه و هم در شبکه‌بندی. تشخیص داده 
شد. رسم هندسه‌ی کانال در این نرم‌افزار در چهار 
مرحله (مختصات‌دهی نقاط رسم خطوطء تخصیص 
سطوح و ایجاد حجم) صورت گرفت. مشخصات کانال 
مربوط به آزمایش 103 (تومیناگا و همکاران) در 
جدول (۱) ارائه شده [4] و در شکل (۲) هندسه‌ی 
کانال ساخته شده در محیط نرم‌افزار ظ0۳ 10511 
نشان داده شده است. پس از این مرحله حجم ساخته 
شده به‌کمک همین نر‌افزار شبکه‌بندی گردید. 


جدول ۱ مشخصات هندسه‌ی کانال موردنظر 
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بهکا رگیری دینامیک سیالات محاسباتی در مدل‌سازی ... 


شکل ۲ هندسه‌ی کانال در محیط نرم‌افزار 0۳ 10171 


مش‌بندی محدوده‌ی حریال. دقت حل مسأله به تعداد 
سلول‌های موجود در شبکه بستگی دارد و اصولا 
شبکه‌های غیریکنواخت. با تعداد سلول‌های مختلف در 
قسمت‌های مختلف بدنه بسته به موقعیت آن, در حبل 
یک مسأله قابلیت بیش‌تری دارند. البته باید درنظر 
داشت که ایجاد چنین شبکه‌ای بستگی به عوامل 
مختلفی از جمله نوع مدل اغتشاشی» شکل ظاهری 
مدل مانند مقطع عرضی نوع مسأله‌ی مورد بررسی و 
مواردی از اين قبیل دارد. رفع خطاهای موجود در حل 
مسأله از طریق تغییر در شبکه و بهینه نمودن آن 
صورت می‌پذیرد و در اکثر موارد تغییر شبکه شامل 
ریزتر نمودن آن در قسمت‌های مورد نیاز است تا 
جایی که نتایج حاصل مستقل از ابعاد سلول‌های موجود 
در شبکه شوند و نتایج کلیدی تغییر نکنند. این امر به 
کمک سعی و خطا و تکرار حل در دفعات مکرر با 
شبکه‌بندی‌های متفاوت حاصل می‌شود. 

در این مطالعه برای شبکه‌بندی محدوده‌ی جریان 
نیز از نرم‌افزار 0۳ 101 بهره گرفته شده است. این 
نرم‌افزار قابلیت ایجاد دو نوع شبکه‌ی چهاروجهی 
(0۵60721ع۲) و شش‌وجهی ([۳۲6:۵۳6۵0۲۵) را دارا 
است و این امکان را در اعتیار کاربر قرار می‌دهد که 
در قسمت‌های مختلف تعداد گره‌ها و المان‌ها را 
متناسب با هندسه تنظیم نماید. از آن‌جایی که مسأله‌ی 
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موردنظر به‌شدت نسبت به شبکه‌بندی محدوده‌ی 
جریان به‌ویژه قسمت سطح آب و دیواره‌های کانال 
حساس است و شبکه‌بندی نزدیک سطح آب در نتایج 
حاصل از تحلیل کاملاً مثر می‌باشد. به‌کمک قابلیت 
این نرم‌افزار در ایجاد این نوع شبکه هندسه‌ی کانال 
شبکه‌بندی گردیده است. 

در این راستا تعداد گره‌های موجود در جهات و 
قسمت‌های مصلف: کاتال فاوت ات و ید از 
چندین مرحله حساسیت‌سنجی نتایج نسبت به 
شبکه‌بندی» شبکه‌ی مناسب انتخاب می‌گردد. این شبکه 
از نوع شش‌وجهی و در قسمت دیواره‌ها و به خصوص 
سطح آزاد متراکم‌تر است که این مسأله از زیاد بودن 
تعداد گره‌ها در این نواحی نشأت می‌گیرد. تعداد گره‌ها 
و المان‌های موجود در کل حجم به‌ترتیب ۱۰۷۹۱۰۰ و 
۸ می‌باشد. برخلاف محدوده‌ی عبور جریان 
آب, شبکه‌بندی در قسمت بالای سطح آب که شامل 
هوا است درشت‌تر می‌باشد و این مسأله به عدم 
دخالت آن در نتایج بازمی‌گردد. این عمل از یک طرف 
می‌تواند موجب پایین آمدن سرعت همگرایی و از 
طرف دیگر به‌دلیل کاهش تعداد المان‌ها موجب افزایش 
سرعت حل شود. در شکل (۳) نتایج مربوط به 
حساسیت‌سنجی مش در شکل گیری سلول‌های جریان 


ثانویه ارائه شده‌است. 
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ب) تعداد المان‌ها -۵۴۸۴۱۸ 


د) تعداد المان‌ها۳۱۶ ۲۰۱۵ 


چ ۲ ۲ 7 7 


۳۱ 


و 6 ۵ ۱۳ 
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ی 


__ ۱ ۷ 
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الف) تعداد المان‌ها- ۵۱۳۷۰۸ 


ج) تعداد المان‌ها ۱۰۴۷۷۰۸ 


شکل ۳ حسایت‌سنجی مش در شکل گیری سلول‌های جریان ثانویه و انتخاب شبکه‌بندی با تعداد ۱۰۶۷۷۰۸ المان 


شرایط مرزی. با توجه به آن‌که تمام سیستم‌های 
فیزیکی در یک فاصله‌ی معین دارای مرز می‌باشند» در 
محاسبات باید شرایط مشخصی در حدود یا مرزهای 
سیستم فیزیکی در نظر گرفته شوند. این شرایط, شرایط 
مرزی نامیده می‌شوند [ ۱۵ ]. اکشر شرایط مرزی در 
کانال‌های روباز شبیه به شرایط مربوط در مجراهای 
بسته هستند. تنها شرط مرزی در سطح آزاد بین این دو 
دسته از جریان‌ها متفاوت است [1۵]. 

از دیگر شرایط مرزی. شرایط مرزی مربوط به 
ورودی و خروجی است. در صورتی‌که جریان فوق 
بحرانی باشد. به هنگام اعمال شرایط مرزی باید عمق 
جریان در ورودی و خروجی مشخص باشد. اما او 
جریان زیر بحرانی باشد. باید عمق جریان در خروجی 
معین و از حل عددی عمق جریان در ورودی محاسبه 
شود و یا این‌که عمق جریان در ورودی معین و از حل 
عددی عمق جریان در خحروجی محاسبه شود [1۵]. 
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شرایط اعمالی به جداره‌ها نیز جزو شرایط مرزی 
محسوب می‌شوند و لازمه‌ی حل معادلات می‌باشند. در 
زیر به شرایط مرزی اعمال شده برای مدل‌سازی عددی 
موضوع مورد نظر در نرم‌افزار 0۳ ۸۱5۷5 پرداخته 
خواهد شد. 


شرایط مرزی در ورودی کانال (0160). با توجه به 
این که جریان در مدل مورد نظر زیربحرانی 
(۳۳-0.46>1) است. لذا توجه به نکات یاد شده در 
بالاه محاسبه‌ی عمق جریان در ورودی به نرم‌افزار 
واگذار شده و از شرط مرزی مقدار مشخص برای 
سرت ره آشبه‌غن نارق زمرق 020 1۷0 و 
5 ۷۷2) در این مقطع استفاده شده است. 
مقدار سرعت در این مقطع بر اساس محاسبات تجربی 
انجام شده توسسط تومیناگاو همککاران [4] 
0.3733 - ۷۷ (سرعت در طول کانال). در نظر 


گرفته شده است. این شرط مرزی دز جریان‌های 
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۳۲ 


تراکم‌ناپذیر مورد استفاده قرار می‌گیرد. در جریان‌های 
تراکم‌پذیر استفاده از این شرط مرزی به نتیجه‌های 
غیرفیزیکی منتهی خواهد شد. زیرا این شرط ورودی 
اجازه می‌دهد که شرایط سکون به هر اندازه‌ای افزایش 
یابد [۸]. 


شرایط مرزی در خروجی کانال (0۷664). در 
خروجی میدان جریان فرض می‌شود که تغییر کمیت‌ها 
(از جمله ارتفاع آب) در راستای جریان قابل اغماض 
است و خطوط جریان به‌صورت مستقیم و موازی 
یکدیگر هستند. بنابراین در این مقطع برای کلیه‌ی 
کمیت‌ها در راستای جریان گرادیان صفر در نظر گرفته 
می‌شود. پس از انتخاب شرط مرزی سرعت در ورودی 
و با توجه به زیربحرانی بودن جریان» شرط مرزی که 
برای خروجی کانال در نظر گرفته شد. عمق جریان 
اشسبت: تفر« قافتا ارو سیک (امتاد 
۵ به این مقطع اعمال گردیده است. همان‌طور 
که فتاه کر ال عفق اب در کانان رایس مستار تانست 
0 0.0905 < ۷ می‌باشد. 

شرط مرزی خروجی فشار احتیاج به‌بیان فشار 
(نسبی) استاتیک در مرز خروجی دارد. مقادیر فشار 
استاتیک معلوم فقط وقتی به‌کار می‌رود که جریان 


زیرصوت باشد. [۸. 


شرایط مرزی جداره‌های کاثال (۷۷۵11). شرط مرزی 
دیوار برای محدود کردن نواحی سیال با جامد به‌کار 
می‌رود. در مسأله‌ی موردنظر جداره‌ها شامل کف و دو 
دهاز کاتال شم باتیله کر ان دارههتا ارات 
پذیرش (۷۷۵۱۱ ٩110‏ ۲۷0) جدار بهره جسته شده و برای 
صحت‌سنجی براساس مدل آزمایشگاهی موجود 
مسطح ([[۷۷۵ 5۳000) در نظر گرفته می‌شوند. 


شرایط مرزی سطح آزاد (6۲609. جریان با سطح 
آزاد به یک حالت جریان چندفازی اطلاق می‌شود که 
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بهکا رگیری دینامیک سیالات محاسباتی در مدل‌سازی ... 


فازها به‌کمک یک سطح مشترک مشخص از هم جدا 
می‌شوند. مثال‌هایی از جریان‌های با سطح آزاد شامل 
جریان کانال‌های روباز. جریان حول 5ظ طلطی 
مسائل خاکریزی مخازن و بسیاری از موارد دیگر است. 
جریان‌های با سطح آزاد با مدل‌های همگن باید در 
مواردی که امکان‌پذیر است استفاده شوند. جریاهای 
پاسطح آزاد با مدل غیرهمگن می‌توانند برای پذیرفتن 
جدایش دو فاز جریان استفاده شوند [116]. 

کاربرد جریان با سطح آزاد در ۳26) به‌طور اساسی 
مشابه حالت چندفازی (همگن یا غیرهمگن) با برصی 
گزینه‌های گسسته‌سازی ویژه برای حفظ تیزی سطح 
تشر کت است. آینها شامل موازد ینس ‌باشت 1161 
۱- یک طرح متفاوت فشرده برای جابه‌جایی افقی 
کسرهای حجمی در معادلات کسر حجم. 
۲- یک طرح انتقال فشرده برای معادلات کسر حجمی 
(اگر مسأله شامل انتقال است). 
۳- رفتار ویژه‌ی ترم‌های گرادیان فشار و جاذبه برای 
اطمینان از این که جریان در سطح مشرک در رفتار 
وی بیافی. مین‌ماند: 

در این نرم‌افزار برای سطح آزاد جریان از شرط 
مرزی 00601718 استفاده می‌شود. نرم‌افزار (۷3 
شرط 0۵060108 را به‌دلیل وجود دو فاز آب و 
هوا و در سطح مشترک این دو نوع سیال در نظر 
می گیرد. 


انتخاب مدل آش شفتگ 
یک مدل آشفتگی عبارتست از یک رویه‌ی محاسباتی 
برای بستن سیستم معادلات جریان متوسط به گونه‌ای 
که کم و بیش بخش وسیعی از مسائل جریان را بتوان 
حل کرد. برای که مان اکن در موارد کلی. 
برنامه‌ی 079 باید دارای کاربرد وسیج» دقیق. ساده و 
از نظر اقتصادی قابل اجرا باشد [12]. همان‌طور که 
اشاره شدء جریان‌های انویه به‌دلیل ناهمگنی آشفتگی 


ایجاد شده و در نهایت القا شده توسط آشفتگی هستند. 
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بنابراین انتخاب نوع مدل آشفتگی به‌طوری که بتواند 
در نمایش سلول‌های این جریان‌ها موثر باشد از 
مهم‌ترین مراحل مدل‌سازی عددی جریان‌های ثانویه در 
کانال‌ها است. 

بنابراین پس از بررسی و ارزیابی چندین مدل 
آشفتگی شامل مدل‌های ۰1-۶ 55۲ 950 و 1.85 
مدل 5560 برای مدل‌سازی سلول‌های جریان ثانویه 
به‌تر از سایر مدل‌ها عمل کرد و نتایج آن در 
صحت‌سنجی (۷۲[7) مدل قابل قبول‌تر از سایر 
مدل‌ها تشخیص داده شد [۱]. از این‌رو در مدل‌سازی 
سلول‌های جریان در کانال مورد نظر از این مدل 
استفاده شده است. در ادامه به بررسی اجمالی این مدل 
آشفتگی خواهیم پرداخت. 

مدل ۹56 (1ع1۱۷۵0 وععقو علامصرع۲ ۵96) 
یکی از زیرمجموعه مدل‌های آشفتگی معادله‌ی 
تش‌های رینولدز ( 16ع1۷00 وعع5 ۳56۲۲۵۱6 
جمنلع۲0۲) محسوب می‌شود. دقت این مدل برای 
جریان‌های پیچیده بسیار بالا است. همگرایی چرخه‌ی 
بیرونی در این مدل نسبت به مدل 6-5 کندتر صورت 
می‌گیرد و علی‌رغم دقت بالای اين مدل. همگرایی در 
آن بسیار آرام صورت می‌پذیرد [16]. 

این :عدل توسط اسیویال: سار کان ق کاسکنی. [17] 
توسعه یافته‌است و از یک رابطه‌ی درجه دوم برای 
همبستگی فشار- کرنش» استفاده می‌کند. برعسی از 
نواقص مدل‌های تنش رینولدز برای شبیه‌سازی 
لایه‌های مرزی, که از معادله‌ی مربوط به ترم ۶ نشأت 
می‌گیرد. سبب توسعه‌ی یک مدل تنش رینولدز شده 
آشت که در ان به جای معادلات مربوط به ترم 6 از ترم 


۱006۵ اه مرن ]۱ سس 
16 ۱۷۵۱۵۲ سوت 


۳۳ 


استفاده می‌شود. این مدل. همانند دو مدل دنکن 
7 برپایه‌ی معادلات انتقال برای همه‌ی مولفه‌های 
مدل فرضیه‌ی ویسکوزیته‌ی ادی را استفاده نمی‌کند اما 
یک معادله برای تنش‌های رینولدز در سیال به‌کار 
می‌گیرد [16]. مدل ۹50 به‌شدت نسبت به شرایط 
شبکه‌بندی در نزدیکی دیواره حساس می‌باشد و برای 
شبیه‌سازی جریان‌های ثانویه مدل مناسبی است. 


نتایج حاصل از مدل‌سازی عددی 


در ادامه نتایج حاصل از مدل‌سازی عددی در محیط 
نرم‌افزار ارائه گردیده است و همان‌طور که در ابتدا 
اشاره شد برای ارزیابی مدل ساخته شده از نتایج مدل 
آزمایشگاهی تومیناگا و همکاران [4] بهره گرفته شده 


است. 


سط حآزاد. شکل (4) وضعیت سطح آزاد را در قالب 
نمودار نمايش می‌دهد. همان‌طور که در این شکل 
ملاحظه می‌شود. در ابنتدای کانال و نزدیک ورودی: 
سطح آب نسبت به مقدار واقعی آن نوسان دارد و با 
حرکت در طول کانال و رسیدن به طول توسعه‌یافتگی 
(تقریبا از فاصله‌ی 1/۵ متری از ورودی), سطح آب به 
مقدار ثابت و صحیح آن نزدیک می‌شود و خطای 
محاسباتی که در ابتدای کانال وجود داشت. از بین 
می‌رود. لذا می‌توان الگوریتم به‌کار رفته در نرم‌افزار را 
برای محاسبه‌ی سطح آب مناسب دانست و به آن 


استناد نمود. 


۲1۲6686 


شکل ۶ نمودار سطح آزاد آب در کانال 
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شکل * سلول‌های جریان ثانویه در آزمایش تومیناگا و همکاران [4] 
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آم)ز۲عجون ]۱ اس 
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شکل ۷ مقایسه‌ی سرعت عمق میانگین بين نتایج مدل‌سازی عددی. آزمایشگاهی و تحلیلی 


سلول‌های جریان ثانویه. شکل‌گیری سلول‌های 
جریان انویه تقریباً از فاصله‌ی ۳ متری از ورودی کانال 
روبه تکامل می‌نهند و در فاصله‌ی ۶ متری از ورودی 
کانال شکل واضحی از این سلول‌ها مشاهده می‌شود 
شکل (۵). تا فاصله‌ی حدود 1/۵ متری از ورودی 
شکل این سلول‌ها متغیر است و ابعاد آن‌ها تغییر می‌کند 
و تقریباً از این فاصله به بعد شکل سلول‌ها ابت 
می‌شود و مانند آنچه که در شکل (۵) نشان داده 
شده‌است» دیده می‌شوند. این مسأله را نیز می‌توان به 
توسعه‌یافتگی جریان پس از طی فاصله‌ی 1/۵ متر از 
ورودی مربوط دانست که مطابق با نتایج آزمایش‌های 
تومیناگا و همکاران است (آن‌ها طول لازم برای 
توسعه‌یافتگی جریان را 7/6 متر دانستند). شکل (") 
نتایج کار آن‌ها را نشان می‌دهد. همان‌طور که در این 
دو شکل دیده می‌شود. نحوه‌ی شکل‌گیری. جهت 
چرخش اندازه و ابعاده محل قرارگیری سلول‌ها و در 
نهایت حرکت به سمت گوشه‌ی آن‌ها مشابه یکدیگر 
است و موید صحیح بودن نتایج حاصل از این بررسی 
عددی است. ابعاد مربوط به دو سلول نزدیک دیواره‌ی 
مایل و محل قرارگیری آن‌ها در فاصله‌ی 1/۵ متری از 
ورودی کانال کاملاً بر نتایج آزمایشگاهی منطبسق 
می‌شود و اين مسأله را نیز می‌توان دلیلی بر صحت 
عملکرد الگوریتم به‌کار رفته در مدل‌سازی سلول‌ها 
#اشگن 

همان‌طور که قبلا اشاره شد حرکنت سلول‌هاق 
جریان‌های ثانویه در کانال‌های ذوزنقه‌ای به‌سمت 
گوشه‌ها قابل توجه است و در نتایج مدل‌سازی عددی 
نیز همان‌طور که در شکل (۵) نشان داده شده‌است. 


سال بیست و پنجم. شماره یک ۱۳۹۲ 


پن پدیده کاملاً مشهود است. اما حرکت روبه گوشه‌ی 
ین سلول‌ها در مدل عددی بیش‌تر از مدل 
آزمایشگاهی است. اینن موضوع را از یک طرف 
می‌توان به کم بودن تعداد نقاط در اندازه‌گیری 
آزمایشگاهی مربوط دانست (۱۰۰ نقطه‌ی اندازه گیری 
لگوریتم به‌کار رفته در آن برای محاسبه‌ی سلول‌های 
جریان. مربوط دانست. اگرچه. بحث درباره‌ی توانایی 
و عدم توانایی نرم‌افزار نیاز به تسلط کامل بر آن دارد؛ 
زیرا به‌عنوان مثال مدل اغتشاشی مورد استفاده ضرایبی 


دارد که هر یک از این ضرایب به نوبه‌ی خحود ار 
ویژه‌ای روی شرایط جریان دارند و کار روی این موارد 
خود بحث دیگری را می‌طلبد. 


سسرعت متوس طگیسری‌شسده‌ی عمقی. برای 
صحت‌سنجی نتایج عددی مربوط به سرعت عمق 
متوسط از نتایج کار آزمایشگاهی تومیناگا و همکاران و 
هم‌چنین نتایج کار تحلیلی شینو و نایت (516) [6] 
استفاده شده است. همان‌طور که در شکل (۷) مشاهده 
می‌شود. نمودار سرعت عمق میانگین محاسبه شده 
براساس سرعت‌های عمقی تست انم فه‌ کمک 
نرم‌افزار تطابق بسیار خوپی با نتایج کار آخهتا خاشتته و 
فقط در نزدیکی گوشه‌ی کانال ( 0.05-> 72 >0.1-) 
مقداری خطا در حدود ۵/ وارد محاسبات شده است. 
به‌نظر می‌رسد این مقدار خطا مربوط به حرکت 
سلول‌های جریان انویه به سمت گوشه‌ها می‌باشد که 
در مدل‌سازی عددی حرکت سلول‌ها به سمت گوشه 


پیش تر از مدل تجربی آن است. 
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۳۹ 


۲ اما 


0۱۱۱02 506۵۲ ۹ 


بهکا رگیری دینامیک سیالات محاسباتی در مدل‌سازی ... 
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شکل ۸ مقایسه‌ی تنش برشی کف کانال بین نتایج مدل‌سازی عددی؛ آزمایشگاهی و تحلیلی 


اساسی در بررسی ارات جریان روی جداره‌های کانال 
است و به دلیل مشکلاتی که در پی دارده از حمله 
مسائلی است که بسیار مورد توجه قرار دارد. از این‌رو 
در این مقاله یکی دیگر از مواردی که کارایی دینامیک 
سیالات محاسباتی و مدل‌سازی عددی به‌کمک آن مورد 
ارزیابی واقع شده است؛ تنش برشی است. 

برای ارزیابی در اين زمینه نیز از همان دو مدل 
آزمایشگاهی و تحلیلی استفاده شده است. در ادامه 
مقایسه‌ای بین نتایج حاصل از مدل‌سازی عددی و مدل 
آزمایشگاهی و مدل تحلیلی شینو و نایت [6] صورت 
گرفته و نتایج حاصل از مقایسه در شکل (۸) نشان داده 
شده اننعتا: مقایسه‌ی صورت گرفته در این شکل. 
حاکی از انطباق نسبتاً خوب نتایج مدل عددی با مدل 
آزمایشگاهی و تحلیلی در ناحیه‌ی کف و انطباق بسیار 
خوب این نتایج در ناحیه‌ی دیواره‌ی کانال می‌باشد. 
البته بیش‌ترین اختلاف بین نتایج عددی و آزمایشگاهی 
مربوط به نواحی نزدیک گوشه‌ی کانال با حداکثر 
خطایی در حدود ۱ بوده‌است و این اخحتلاف نیز به 
همان عوامل توضیح داده شده در بخش قبل مربوط 


ین کرکف: 


نشریه مهندسی عمرآن فردوسی 


نتیجه گیری 

در حالت کلی نتایج تهوست امه از دل‌سازی توق 
سلول‌های جریان در کانال را به‌صورت زیر می‌توان 
ارائه نمود: 

روش عددی و الگوریتم مورد استفاده برای 
مدل‌سازی سلول‌های جریان انویه در این نرم‌افزار 
مناسب بوده‌است و تطابق بسیار خوبی بین نتایج مدل 
عددی و آزمایشگاهی در نمایش سلول‌های جریان 
وجود دارد. لذا درکل می‌توان دینامیک سیالات 
محاسباتی را توانمند در مدل‌سازی این سلئول‌ها 
واتنیت: 
۱- مواردی از قبیل تعداد المان‌ها. نوع مدل اغتشاشی. 
شرایط مرزی و ... در مدل‌سازی سلول‌های جریان 
بسیار موثر است. 
۲- مدل اغتشاشی به‌کار گرفته شده (تنش رینولدز 
راشای تاک بارخ ماسی ی تیا 
۳- نمودار سرعت عمق متوسط تطابق بسیار خوبی با 
نتایج آزمایشگاهی و تحلیلی دارد. 
6- مقایسه‌ی نتایج حاصل از محاسبه‌ی عددی تنش 
برشی مرزی با نتایج آزمایشگاهی و تحلیلی حاکی از 
انطباق خوب نتایج به‌ویژه در ناحیه‌ی دیواره‌ی جانبی 


انشدتا: 
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۳ ورستیگ. هنگ کارل.. مالالاسکرا؛ ویراتونگ «مقدمه‌ای بر دینامیک سیالات محاسباتی»» ترجمه‌ی شجاعی‌فرد. محمدحسن. 


نورپورهشترودی» علیرضاء انتشارات دانشگاه علم و صنعت. چاپ دوم. تهران» (۱۳۸۰). 
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۲۸ بهکا رگیری دینامیک سیالات محاسبانی در مدل‌سازی ... 


م28 ۷۵۱0 ۱۱6۵-۸۷۵۵۵60 مصنلع۱۷]۵0 ربعم رع120 200 ,۱۷۲ رصن ,۷۷ .1 راطعتصک .14 
۰ ,11۳۱9۱۱66۳۱۴۱9 30۲۵۱6 ]۵0 دیول رما ته0ج0۵ع5 طزمه کلعصصف اهلذ۲۳۵۵۵2۵ صز تهعدو 
2007(۰) ,)0۸۹5 ر1 .0 ر133 
۵ گودرزی. محسن و عظیمیان. احمدرضاء «مدل‌سازی عددی جریان در کانال روباز به‌روش ردیابی سطح آزاد», مجله‌ی دانشکده 
فنی دانشگاه تبریزه ص 0۷-۶۵ (۱۳۹۱). 
2009(۰) را ۳۲0۵6۰ ۸۳۲5۷5 ,13 ,۸۳۲۷5 ,712.0 ,افتاصعطظ تعوتا 0۳2 ۷5 5]ظظ .16 
۶ صمتاحاعنی صتحتاومتناوومم معط فصتامل۷۵ ریظ.] رتکاعاهت) 200 بو وتقه9 وتا ر6ا212ع(5 .17 
-245 .00 ,277 ۷۵ رک۷6۵۳۱۵۳۱( یا .ل موه عصماوری آق‌نصمصر اصمهصا چم بععصملتا تنم 


272, )1991(. 
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